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We present a general method for systematically finding the families of exact ground states of 
classical Heisenberg models with the following properties: (i) they have the translational symmetries 
of infinite crystals; and (ii) they contain only bilinear terms, i.e. Si-Sj. With very limited exceptions, 
the method is applicable to all such models with one spin per unit cell, and models that have many 
spins per unit cell and satisfy certain conditions. Crystal dimensionality, range of coupling, etc. 
do not affect the validity. We also prove theorems that guarantee the existence of, in certain 
circumstances, a spiral state within the ground state manifold. Our formalism is predictive for 
disorder, and it might be generalized to higher-order terms (i.e. (Si • Sj) r ), anisotropic couplings, 
and arbitrary dimensional spins. 
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Exchange interactions dominate magnetic orders in 
solids. Despite the quantum nature, many studies on 
magnetic properties were set out classically. The rea- 
son, besides simplicity, is that classical models are indeed 
good approximations to quantum ones in certain cases, 
two of which are of large spins and of spin— ~'s. The 
former has a footage in large-spin expansion [1], whereas 
the latter relies on the Hartree-level equivalence between 
spin— | 's and classical spins. Large localized spins ap- 
pear quite commonly in rare earth materials. To achieve 
quantum anomalous Hall effect (QAHE) on them, a non- 
coplanar classical spin order is required so that the T- 
reversal symmetry is broken On the other hand, in 
quest for quantum spin liquids (SL), investigations are 
often carried out as a quantum fluctuation analysis that 
succeeds classical calculations (see e.g. [3]). In these and 
many such ongoing fields, classical calculations play a 
fundamental and crucial role. 

It is notoriously nontrivial to deduce the magnetic or- 
ders from a formally simple exchange interaction Hamil- 
tonian. While numerous profound techniques have been 
developed in the quantum scenario (e.g. mean field the- 
ory large-spin expansion H, spin wave theory 
high temperature expansion |6[, transfer matrix meth- 
ods 0, [sf? and renormalization group theory [9]), less 
attention has been paid to classical models. To name a 
few, one prevalent idea is local minimization: instances 
are termwise [10] or cluster [ll[ minimization, as well as 
rewriting the energy in special local forms (see e.g. [12|). 
Another pop ular device is the classical spin wave theory 
(see e.g. Hall!), 

which usually assumes a spiral or he- 
lical ansatz (see Fig. [1]) within which the energy is min- 
imized. Furthermore, the Luttinger-Tisza (LT) meth- 
ods first weaken the normalization constraints on spins 
and subsequently check the strong constraints (see e.g. 

There are also methods founded on some classical- 
quantum correspondences (see e.g. [HI). These, how- 



ever, are not enough to give us a comprehensive under- 
standing of classical spin orders, and a more satisfactory 
method is demanded. 

In this Letter, we present a general method for finding 
the family of exact ground states of a classical exchange 
interaction Hamiltonian with the following properties: (i) 
the Hamiltonian has the translational symmetry of an 
infinite crystal; and (ii) the Hamiltonian contains only 
Heisenberg terms, i.e. Si • Sj with the spins being real 
unit 3- vectors (these will be referred to as the two "basic 
assumptions"). We emphasize on the generality of the 
Hamiltonian: We did not make any assumption on the 
crystal dimensionality, the particular Bravais lattice, the 
internal structure of unit cells, or the range of exchange 
interactions. Thus it can be square, Kagome, face cen- 
ter cubic, or diamond structures with interactions among 
nearest-neighbors, second nearest-neighbors, and so on. 
We shall prove that, if the system has only one site (spin) 
in each unit cell, or if a system has many spins per cell 
but satisfies certain explicit conditions, then among the 
ground states there must exist a spiral state (see Fig. [1]). 
We shall also show that, with fairly limited exceptions, 
our method can find all the ground states of any such 
system. Our "ground states" are not computed within 
any ansatz; namely, they are rigorous. 

This Letter is organized as follows. We first focus on 
systems with one spin per unit cell; we shall prove the spi- 
ral theorem and give the guidelines for finding the whole 
ground state manifold (GSM). Then, with a little addi- 
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FIG. 1: Schematic (a) spiral and (b) helical states. 
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tional complication, we present a similar set of results 
for systems with multi spins per cell. Finally, we com- 
ment on the cases we cannot easily handle, and discuss 
a few other possible generalizations of the framework. It 
may be helpful to bear in mind that our method consists 
of only three essential steps: changing variables, figur- 
ing out the new restrictions, and determining the ground 
states. 

One spin per unit cell. — A general classical exchange 
Hamiltonian in this case takes the form 

Hd = ^ ' Sfi+ffj , (1) 

n j 

where real unit 3- vector Sft is the classical spin at lattice 
point n; ffj connects the two lattice points involved in an 
exchange term with parameter J 3 \ and the ^ is due to 
double count. In this Letter we will keep the convention 
of positive signs before coupling constants. This problem 
is completely classical, but for convenience, we reformu- 
late the energy expression with the help of a Hilbert space 
7-L that has an orthonormal basis {|n)}. Suppose there 
are N spins. Viewing Sft as column vectors, we define the 
"ket state" and "bra state" corresponding to a particular 
spin configuration as 

w = 4=£3i®ifl>, w = iE4®H ( 2 ) 

n n 

with f representing the Hermitian conjugate of a matrix. 
We also define the effective Hilbert space Hamiltonian 
corresponding to the classical H c i as 

H = ^EE^i^+^i- ( 3 ) 

n j 

Then one figures that the quantity 

e={mW) (4) 

recovers correctly the classical energy per spin. The 
bonus for defining H is that in as an operator, H 
has an unambiguously defined set of eigenstates |fc) = 
^= Y^fi elk ' n \n) and eigenvalues eg; we refer to the latter 
as the spectrum of H. One can also show from the ba- 
sic assumptions the property e_g = eg. After change of 
basis, one finds Eq. (j4]) can be written in terms of eg as 

k 

with the new variables 5g being the Fourier conjugates: 
%=7f E^e-^. (6) 

ft 

The above procedure formally looks as if we were defining 
quasi-particle excitations by Eq. (|6j) so as to diagonalize 



the Hamiltonian, but 5g or Sfi are by no means annihila- 
tion operators. Now that we have a "diagonalized" form 
Eq. (|5j), we need to know what values {£g} can take. 
Noticing that the set of allowed {Sfi} is 

V = {{Sn}\SfteR 3 , 4^ = 1}, (7) 

one can prove that a choice of {5g} is allowed if and only 
if it is contained in 

1 = e C 3 , S_ n = St, SlS Up = NS dp ,W}.(S) 

k 

V and X are the domain and the image of the map Eq. 
(|6|), respectively. The p in Eq. (|HJ) is assumed to be lo- 
cated at, and the summations of k in Eqs. ([5]) (|Hj) are 
over, the reciprocal lattice points in the Brillouin zone 
(BZ). In addition, in Eq. (|SJ), St is assumed to be peri- 
odic. Now we are ready to prove the following theorem. 

Theorem 1 (Spiral ground state theorem (one-spin 
cells)). For any classical Heisenberg model that satisfies 
our two basic assumptions, and that has only one spin 
per unit cell, among the ground states there must exist 
a spiral state. 

Proof From the definition of X we know J^g jjStS^ = 1; 

we also know jjStS^ > 0. Thus if S^ only at 
some global minima in the spectrum, then the result- 
ing state must be a rigorous ground state. Suppose £?o is 
a global minimum; since e_g = eg, so is — fco- One can 

verify that the specification S ± ^ q = ^^(l,±i,0) T for 
inequivalent ±ko (i.e., not connected by reciprocal trans- 
lational vectors), or v // ^V(l,0,0) T for equivalent ±fco? is 
contained in X and is hence allowed. More explicitly, 
Sn = (cos(£ • n),-sin(fe ■ ™),0) T or (e^° ,f \ 0, 0) T , for 
the two cases respectively. These are spiral states. □ 

From the proof one concludes that a sufficient and nec- 
essary condition for an allowed state to be a ground state 
is that its Sg vanish at all not-the- lowest points, so as to 
have the same energy as the spiral state given by Theo- 
rem [TJ This is the principle for finding the entire GSM: 
For any given Hamiltonian, we compute and examine the 
spectrum, write down equations according to X for 5g at 
the minima, and determine the GSM. Moreover, since 
in this one-spin case, the distribution of spectral min- 
ima completely determines the GSM, it can be used for 
classification of all models. Be aware, however, that the 
distribution may be so special that the equations imposed 
by X are modified due to the periodicity of 5g (e.g. Figs. 

|5(e)- (h)) or dzfei being equivalent (e.g. Figs. [5(a)- (d)). 

We shall now classify two-dimensional (2D) systems in 
the afore- mentioned way and present their GSMs (l6| . 
We are free (for any lattice) to assign integer numbers to 
n Xl n y so that the BZ is (— 7r, tt] x (— 7r, tt] (we will use this 
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FIG. 2: For 2D systems with one spin per cell, the special 
distributions of spectral global minima, represented by the 
dots, in the BZ (— 7r,7r] x (— 7t,7t] (see text). The larger and 
the smaller dots are respectively defined as zb/co and zbfci. 
This list is complete for the case of one pair of minima ((a) - 
(h)) but partial for the case of two pairs ((i) - (r)); definitely 
there can be more. In (o) - (r), ±ki are generically located. 



convention throughout the Letter). Special locations of 
the minima are partially enumerated in Fig. [2j the GSM 
of them and of some general cases are summarized in 
Table H The expressions in Table U literally correspond to 
spiral phase, ferromagnetic (F) /antiferromagnetic (AF) 
phases, frustrated phases, spiral phase, frustrated phases, 
and helical phases (possibly with alternating helicities). 
Although the cases enumerated are obviously incomplete, 
one realizes that the GSM of any given Hamiltonian can 
be similarly worked out extremely easily, unless there are 
infinite number of spectral minima. Our formalism also 
makes phase transitions manifest, as there is usually an 
abrupt change in the number or locations of the lowest 
points when a transition occurs. 

Let us apply our method to the J\ — J2 models 
on square and triangular lattices, where Ji, J 2 cor- 
respond to nearest- and next nearest-neighbors respec- 
tively (recall our positive sign convention). The spec- 
trum for square lattice is = Ji(cosk x + cosk y ) + 
2 J2 cos k x cos k y . The regimes with J\ < 0, J2 < 0.5|Ji|; 
J\ > 0, J 2 < 0.5Ji; and J 2 > 0.5 respectively fall un- 
der Figs. [2|(a) (d)(1), and hence have F, AF, and frus- 
trated two-sublattice AF orders. The spectrum for tri- 
angular lattice (for some particular translational axes) 
is e£ = Ji[cosk x + cosky + cos(k x — k y )\ + J2[cos(k x + 
k y ) + cos(2k x — k y ) + cos(2k y — k x )\. The regime J\ < 
0, J2 < ^\Ji\ is given by Fig. [2fa) and is ferromagnetic; 
for J\ > 0, J 2 < there is a single pair of minima 
=b(^7r, — §tt), implying an angle |7r between neighbor- 
ing spins; the remained cases have three pairs of lowest 
points, which are not included in Table U but are solvable. 

Multi spins per cell. — For a general system with N 
unit cells that each has m spins, the Hamiltonian is of 
the form 



TABLE I: The GSMs corresponding to different distributions 
of the spectral global minima, for 2D systems with one spin 
per cell. 7 £ R, p £ are arbitrary and preassigned. 



Spectral global minima 



GSM 



cos(/co • ft) 
p I — sin(/co • ft) 


Sn = p\ 

V / 

y/2 cos (ko • ft + j ) cos 7^ 

Sn = P I v / 2sin(fco • ft + f ) sin 7 


cos (ki • ft) 

Generic ±fc , ±ki Sn = p \ - s'm(k 7 • ft) I , i = or 1 



Generic ±ko 
Figs. m*Md) 
Figs. ESeJ-EZth) 



Figs. ETi)^n) 
Figs. EtoJ-Etr) 



'e^ -"cos7' 
Sn = p I e^ 1 ^sin 7 

^cos(fci • ft) sin 7^ 
Sn — P I sin(fci • ft) sin 7 

e^°"*cos7 



where a^, bj = 1, 2, . . . , m label the particular internal 
sites involved in a coupling term. Using an orthonormal 
basis of a Hilbert space we define 



(10) 



In the basis \ka) = ^= J2n e%k ' n \n°)i Eq. (flQ|) can be 
written as H = a b I® |fea) (ft(fe)) 6 (fc6|, where matri- 
ces ft(fc) have well-defined orthonormal eigenvectors £j* a 
and eigenvalues eg a . In this multi-spin scenario, one ob- 
serves a band-like structure, with a being the band label. 
Some universal properties can be proved: h(—k) = h(k)* , 
e -ka = e ka' and f-E, a = £| a (convention). Following a 
similar procedure we find the per unit cell energy: 



(11) 



k,a 



with the new variables, in an abbreviate form, defined as 

(12) 



and 



(9) 



where =5^ ^ — (S'nl, • • • , ^nm)? ^ £ — ('S'jfeiJ • • • j ^kmJ" 
^k fel' • • ■ >fjfem)" The ima g e ° f the ma P Ec i« C2) is 

also straightforwardly computable: 

I = {{%J|^ a GC 3 , s_~ ka = st a , 

E( E ^uA + ^ = NS ^ v «-p}' ( 13 ) 
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where the a on the second line are matrix indices. Now, 
we shall present a similar spiral theorem for multi-site 
unit cells but with an extra assumption, as below. 

Theorem 2 (Spiral ground state theorem (multi-spin 
cells)). For any classical Heisenberg model with m spins 
per unit cell that satisfies the two basic assumptions, if 
there exist points ±fco and a band ao where e^ a assumes 
the global minimum and such that £,r can be chosen 

to take the form f , t = -j=(e ±u \. . . ,e ±it -) T , then 

^±koao \/rn\ / ' 

among the ground states there exists a spiral state. 

The proof is analogous. The spiral state to be found 
is Sfta = (cos(4 • ft - t a ),-sin(fc ■ n - t a ),6) or 
(e*( feo ' n-t °), 0, 0) for inequivalent or equivalent ±ko re- 
spectively. We further present two corollaries. 

Corollary 1. For a system with two sites per unit cell, 
if the 2x2 matrices h(k) contain no Pauli a z component 
for all fc, then there exists a spiral ground state. 

Corollary 2. For a system with m spins per unit cell, 
if for every /c, (ft(fc)) only depends on (a — b) mod m, 
then there exists a spiral ground state. 

A classification in the multi-spin case should take into 
account the £g . Whenever a spiral ground state is guar- 
anteed, the GSM can be determined by focusing only 
on the spectral minima. Unless there are infinite num- 
ber of them, the calculation will be rather simple. It is 
fortunate that some models do fulfill the prerequisite on 
£g . For instance, Corollary [T] is applicable to the hon- 
eycomb J\ — J2 — J3 model, which has been studied by 
classical spin waves [H, [Hj ; our method can determine 
the ground state energy and, in regimes with finite num- 
ber of minima, the whole GSM. Another instance is the 
checkerboard J\ — J 2 model, where Ji, J 2 correspond to 
straight and diagonal edges respectively (see e.g. [17j), 
and whose h(k) has no o z component along \fc x \ = \k y \. 
For < 1, a single minimum appears at k = 0, imply- 
ing a F ( Ji < 0) or AF (Ji > 0) phase. For ^ > 1, the 
minima cover the entire BZ boundary. 

Discussion. — We cannot easily handle the following 
cases: the case of an infinite number of spectral min- 
ima, and the multi-spin case in which the £g condition 
in Theorem [2] is not satisfied. Nevertheless, one can still 
gain some knowledge within our framework. For exam- 
ple, ideally QAHE desires a purely classical Hamiltonian 
with no coplanar ground state, and our results imply such 
a system must go beyond the spiral theorems. Further- 
more, since the dimension of GSM is not greater than 
the number of independent variables (of the same order 
as the number of nonzero ), a disordered system with 
an infinite dimensional GSM must have infinite number 
of spectral minima or violate the £g condition. 



There exist a few possible generalizations of the frame- 
work. To accommodate Hamiltonians with exclusive 
higher-order terms, i.e. (Si • Sj) r one may redefine 
1^) = ^En,afe)® r ® To include anisotropic 

couplings of the form S^J^S 1 -, one may replace the Jjl 
in Eq. ([TQ|) with a matrix Jj. In addition, the spins 
can be generalized to arbitrary dimension, adapting the 
formalism to XY and Ising models, etc. 

Upon completion of this work, we noticed the similar 
conclusions about spiral states given by the LT methods 
(see e.g. [14]), but we have derived them in a differ- 
ent and more general framework. In this framework, it 
is rather simple to determine the entire GSM, classify all 
models, predict phase transitions, and draw hints for non- 
coplanarity and disorder. Our philosophy of keeping the 
full information, as opposed to the LT methods, seems so 
natural that it should have been discovered decades ago. 
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In this Supplemental Material we show explicitly how to determine the entire ground state manifold (GSM) of a 
system within the framework we have developed. We consider two exemplifying cases for two-dimensional systems 
with one spin per cell: (i) the case in which there is a single pair of generically located minima ±fco m the spectrum; 
and (ii) the case in which there are two specially located minima ko = (0,0), k\ = (tt, 0) in the spectrum. 

Consider case (i), which corresponds to the first item in Table I of the main paper. In this case, any ground state 
can have nonzero only at ±fco- Hence, we set 

4=Vf(|)- s -«i=\/f(*:)- (i) 

and zero otherwise. We must check if such a specification is contained in 

Z = {{%)!% e C 3 , = % ]T SlS Up = NS^VP}. (2) 

k 

The first two conditions are fulfilled. The third is effective when p = and ±2/co (up to a reciprocal translational 
vector). The cases of p = 0, 2/cq respectively require 



XV 



•' 2 ' 'y| 2 + |*| 2 = l, * 2 +^ + z 2 = 0. (3) 



The case of p = — 2ko gives a condition conjugate to that of p = 2ko. The set of solutions to Eq. (j3]) would give the 
GSM. Written in terms of real and imaginary parts, namely x = x r + ix" , y = y f + iy" , z = z' + iz" , Eq. ([3]) reduces 
to 

||(^,^z0|| = \\{x\y\z")\\ = -L, (x',y',z') JL {x\y\z"\ (4) 

where we have viewed (x' ', y', z'\ (x", y" , z") as 3-component vectors. Now, apparently the set of solutions to Eq. ((4]) 
is, up to an arbitrary rotation, 

(x',y',z') = (-j=,0,0), (x",y",z") = (0,^=,0), (5) 

whence 

4^ = — \ I • (6) 




The GSM is therefore 



cos(/co • n) 

Sft = p\- sin(fe -ft)], (7) 







where p G SO3 is an arbitrary rotation. 
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Next, consider case (ii), which falls under the fifth item in Table I of the main paper. The periodicity and the 
definition of X together imply that , are real. We thus set 



N I x \ * [W /u 



where x, y, z, u, v, w G R. The p = (0,0), (tt, 0) cases of StS^^ = iV<5^ give two independent constraints: 

||(x,^,z)|| 2 + ||(^, w)|| 2 = 2, (z >2/ ,*) • (u,v,w) = 0. (9) 
The solutions are, up to an arbitrary rotation, 

(x, y, z) = v / 2(cos 7 , 0, 0), (u, v, w) = y/2(0, sin 7,0), (10) 
where 7 G M is another arbitrary parameter. Consequently, 

(cos 7\ _ 
1, § ni =>/N[sm~ ). (11) 

Therefore, the GSM is 

cos 7 

S fl = />|(-l) n jsin 7 |, (12) 



where p G 503, 7 G 1R are arbitrary. 



